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Objectives 


The primary aim of this unit is to introduce and to define precisely the 
concept of a definite integral. 


After working through this unit you should be able to: 


(i) explain the meaning of the term definite integral: 
(ii) evaluate the definite integral of functions of the form 


xH— an? + bx? + ex" (x € R) 


where a, b and c are real numbers, and p, g and r are real numbers 
none of which is equal to —1: 
(iii) estimate the value of the definite integral for a tabulated function: 
(iv) determine the upper and lower error bounds of this estimate, 
including the effects of errors in the tabulated values; 
(v) estimate the value of the definite integral of a function using the 
trapezoidal rule; 
(vi) estimate the number of steps (and the width of interval) required 
to obtain a given accuracy using this rule; 
(vii) estimate the number of decimal places required in the tabulated 
values to achieve this accuracy: 
(viii) calculate the average of a function over a given interval; 
(ix) calculate the volume of revolution of a given curve; 
(x) calculate the distance travelled by a body in a given time, given the 
velocity as a function of time; 
(xi) translate the traditional Leibniz notation for a definite integral 
into the functional notation used in this course, and vice versa. 


N.B. 

Before working through this correspondence text, make sure you have 
read the general introduction to the mathematics course in the Study 
Guide, as this explains the philosophy underlying the whole course. You 
should also be familiar with the section which explains how a text is con- 
structed and the meanings attached to the stars and other symbols in the 
margin, as this will help you to find your way through the text. 
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Bounded Estimates 
for the Area 

under a Curve 
9.1.1; 9.1.2 


The Definite Integral 
for Combinations 

of Functions 

9.2.2 


The Definite Integral 
9.2.1 


Applications 
Volume of a Solid 


Equation of a of Revolution, 

! Straight Line Averages, Velocity 
1 R.B.3 and Distance 
beweeoccccus 9.3.1; 9.3.2; 9.3.3 
peemonacssrs oa] Approximation Methods 

! Approximating Trapezoidal Rule 

1 Polynomials Simpson’s Rule 

1 Unit 4 | 9.4.1; 9.4.2 


Glossary 


Terms which are defined in this glossary are printed in CAPITALS. 


ABSOLUTE ERROR 


ABSOLUTE ERROR 
BOUND 


AREA OF A 
RECTILINEAR 
FIGURE 


AREA BENEATH THE 
GRAPH OF f 
BETWEEN 

a AND b 


AVERAGE 


BOUND 


DEFINITE INTEGRAL 


END-POINTS OF A 
DEFINITE INTEGRAL 


INTEGRATE 
ORDINATE 
RECTILINEAR 


FIGURE 
SIMPSON’S RULE 


SUB-INTERVAL 


TRAPEZIUM 


vi 


The ABSOLUTE ERROR in a measurement x is the 
difference x — ¥ between the measured number 
and the exact number Y. 

The ABSOLUTE ERROR BOUND in a measurement is 
the maximum possible value of the magnitude of 
the ABSOLUTE ERROR. 

The area of a rectangle is base x height. 

The area of a right-angled triangle is 

4x base x height. 

The area of a RECTILINEAR FIGURE can be deduced 
from these definitions. 

The AREA BENEATH THE GRAPH OF f BETWEEN 
AND b is the limit of the sequence defined by 
S, = All f(a)| + [f(a + hy t+ +1 f(a t+ {a — LAH. 


| Took Seat eer 
where h = ——. if this limit exists. 
n 


The AVERAGE of f(x) over the interval [a, 6] is H. 


where 
1 b 
b-a f f 


The number B is an UPPER BOUND ofa set of numbers 
Sifs < B for each element se S. 

The number d is a LOWER BOUND of a set of numbers 
Sifs > 6 for each element se S. 

A set is BOUNDED if it has upper and lower bounds. 
The DEFINITE INTEGRAL of the function f between 
a and 6 is the limit of the sequence defined by 
S, = Alfa) + fla +h) +--+ + fla + {n -— 1h] 


| a ne re 
where h = = if this limit exists. 
it 


H= 


The END-POINTS OF THE DEFINITE INTEGRAL of the 
function f over the interval (a, 6] are a and b. The 
UPPER END-POINT is b, and the LOWER END-POINT 
is a. 

To INTEGRATE means to evaluate a DEFINITE 
INTEGRAL. 


The orDINATE of a point in the Cartesian plane is 
the y-co-ordinate of the point. 

A RECTILINEAR FIGURE % a figure bounded by 
straight lines. 

SIMPSON’S RULE is an approximation to the DEFINITE 


INTEGRAL OF { BETWEEN @ AND 0 given by the 
formula 


e h 
SL = 00 + 4Y1 + yn to + AYa ay + Ia) 


where the interval [a, 6] is divided into n (n even) 
SUB-INTERVALS and yo, ),,....3, are the ORDIN- 
ATES at the ends of these sub-intervals. 

A SUB-INTERVAL of (a, 4] is an interval [c, d] where 
age<dc¢b. 

A TRAPEZIUM is a figure with four straight sides, of 
which two are parallel. 


Page 
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TRAPEZOIDAL RULE 


VOLUME OF 
REVOLUTION 


WIDTH OF AN 
INTERVAL 


The TRAPEZOIDAL RULE is an approximation to the 
DEFINITE INTEGRAL OF { BETWEEN a AND b given by 
the formula 


2 h 
[f= 300+ a + t+ es +90) 
pra 


where the interval [a, 6] is divided into » suB- 
INTERVALS and yo, ),,...,), are the ORDINATES 
at the ends of these sub-intervals. 

A VOLUME OF REVOLUTION is the volume of a solid 
bounded by a surface, the surface being generated 
by rotating the graph of a function fin an interval 
(a, 6] about the x-axis. 

The WIDTH OF THE INTERVAL [a, d] is b — a. 


vil 


Page 
50 


40 
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Notation 


The symbols are presented in the order in which they appear in the text. 


e, The absolute error in a measurement x. 

& The absolute error bound. 

a, The sum of the areas of n rectangles, which gives an 
estimate of a required area which is less than the required 
area. 

a The sequence a,,a2,43,... 


The sum of the areas of n rectangles, which gives an 
estimate of a required area which is greater than the 
required area. 


A The sequence A,, Az, A3,... 
A The common limit of g and A. 
h The width of a sub-interval of [a,b] when [a,b] is 
divided into n equal sub-intervals. That is, 
h= ee (neZ*) 
n 

S, The sum h[ f(a) + f(a +h) +--- + f(a t+ {n — L}h)). 

b 
f f The definite integral of f in [a, b). 

a 
Bibliography 


Two books which introduce integration before differentiation are men- 
tioned below. The approach to the definite integral in both is more 
rigorous than the approach in our course, and would be useful reading 
for those who wish to extend themselves at this stage. However, do not 
be surprised if you find them a little difficult. 


T. M. Apostol, Calculus Vol. J (Blaisdell 1967). See Part | of the Introduc- 
tion, and also Section 1. 


B. Hunt, Calculus and Linear Algebra (W.H. Freeman 1967). See Chapters 
8 and 9. 
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9.0 INTRODUCTION 


The definite integral is a mathematical concept which arose from the 
problem of determining area. This problem links together two major 
ideas from previous units — inequalities and limiting processes —to 
produce a mathematical formulation of the meaning of area. From this 
we obtain the concept of the definite integral, which we find has many 
other applications. 


If you have met the topic of integration before, as well as the topic of 
differentiation which appears later in our course. you may be wonder- 
ing why we introduce the definite integral first. It is perfectly possible to 
introduce these topics in either order, and each approach can be made 
logically consistent. In fact, we shall introduce integration and differentia- 
tion each in its own right, to emphasize that each is an important mathe- 
matical concept. In a later unit, we shall establish the connection between 
them in the Fundamental Theorem of Calculus. We begin with the definite 
integral because historically the basic ideas of this were known long before 
those of differentiation, although the rigorous treatment of the integral 
came at a much later stage. 


If you have some familiarity with the calculus you may also wonder 
why here, and in the later units on integration and differentiation, we 
use a notation which is different from the classical notation with which 
you are familiar, and which is used in most textbooks. (In fact, you will 
find that there are a number of recent textbooks which use the same 
notation as ours.) There are several reasons for this change. The notation 
we use is consistent with our basic approach to mathematics through the 
concept of a function. Also, we wish to avoid some of the conceptual 
difficulties experienced by many students, which may in part be due to 
the traditional notation. If you are new to the subject there should be 
no notational difficulty. If you have already studied calculus, we suggest 
that you make a fresh start: the new notation should enable you to 
concentrate on the principles (as opposed to the techniques), because of 
its very novelty. Once you have mastered the basic principles, there is 
no objection to your using the traditional notation : the conversion from 
one notation to the other should not prove difficult. 


The history of the problem of finding area is an interesting one. In the 
early years of ancient Babylon it was believed that the area of a plane 
figure depended on its perimeter. (If you were an Estate Agent in Babylon, 
how would you buy and sell your land to take as much advantage of this 
belief as possible?) 
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However, the correct methods for finding areas of rectangles and triangles 
were known before 2200 z.c. The next step, that of finding areas of plane 
figures with specific curved boundaries, such as a parabolic arc, was 
apparently not taken until the time of Archimedes (287-212 B.C.). His 
mathematical method is essentially the same as the one we develop in 
this unit. One of his elementary checks was to cut out the appropriate 
shape in a material of uniform density, and compare its weight with that 
ofa shape of the same material and of known area. Archimedes realized 
that the problem of the determination of volumes bounded by curved 
surfaces was similar to the problem of determination of area, and that 
both could be approached by using a process involving closer and closer 
approximations. Much later, in the seventeenth century, Newton, amongst 
others, formalized integration and established its link with differentiation. 


Newton’s work aroused enormous interest, and the names of many 
mathematicians of renown appear in the history of the calculus. By the 
time Riemann published his definition of the definite integral in 1854, 
mathematicians had long realized that they had in their handsan extremely 
powerful tool. 


In this unit we develop the theory of integration along the following lines. 
We first consider ways of estimating area, and we note that elementary 
methods of calculating areas precisely can only cope with regions bounded 
by straight lines. Other regions require a different approach : we approxi- 
mate to the region by successive sets of rectangles (or other suitable 
shapes of known area), and thus obtain a sequence of what we intuitively 
consider to be better and better approximations to the area. The limit 
of this sequence of approximations is called a definite integral. This 
concept can be applied to a far wider range of problems than the calcula- 
tion of areas. For example, integration is used in determining volumes, 
and rate of flow ofa liquid through a pipe, as you will see in the television 
programme associated with this unit. 


We shall see later that, even though we can represent an area by a definite 
integral, there are cases when we cannot evaluate it in practice. In these 
circumstances we are forced to return to estimating the area ; we conclude 
the unit with a short discussion of two methods of doing this. 


You may have gathered from what we have been saying that we are now 
going to tackle this problem with the guns we have been priming in 
earlier units. We shall, in particular, be using the ideas of function, 
accuracy and errors, and limit of an infinite sequence. The next few units 
should help to show the importance of some of the ideas in the earlier 
units. 


i 


FM 9.0 


Archimedes (Mansell) 


9.1 THE USE OF BOUNDED ESTIMATES IN 
FINDING AREA 


9.1.0 What is Area? 


Intuitively we all know what we mean by area, just as we feel we know 
the meanings of length, time, speed and volume. We are also aware 
that we use each of these words in two senses: sometimes to mean a 
physical quantity, and sometimes to mean a measurement of that quantity. 
For example, the word “‘area’” means a vacant piece of level ground; 
but to avoid using a clumsy phrase, we say “‘the area is five acres”’, where 
the word “‘area’’ stands for “the measurement of the area”. The sense 
intended is usually obvious from the context. Here, we are interested in 
the measurement of area. We can leave the philosophers to worry about 
the fundamental concept. As mathematicians, we use our intuition to 
guide us to a precise mathematical formulation of area; we then develop 
and generalize the mathematical concept we have defined. You have 
already experienced such a development in Unit 7, Sequences and Limits I, 
where we discussed a possible mathematical formulation of the intuitive 
concept of velocity at an instant. 

Using our intuition, let us define* the area of a rectangle as the product 
of the lengths of two adjacent sides. With this as our starting point, we 
shall proceed to define areas of other figures of increasing complexity, 
always assuring ourselves that our definitions accord with our intuitive 
expectations. 

Let us start with rectilinear figures, which are figures bounded by straight 
lines. By fitting together two identical right-angled triangles to form a 
rectangle, we find that the area of such a triangle is half the product of 
the length of a base and the corresponding height, which we shall write 
as 4 x base x height. 


Knowing this, we can find the area of any rectilinear figure by regarding 
the figure as being built up of rectangles and right-angled triangles. 


Exercise | 


* We define area such that it has dimensions L? (see Unit 3, Operations and Morphisms, 
page 36). We shall state the units of area when they are known: otherwise we shall only 
give the magnitude of the area. 
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Discussion 


Definition 1 


Definition 2 


Exercise 1 
(2 minutes) 


A trapezium is a quadrilateral which has two parallel sides. 
Show that the area of the particular trapezium shown is 3ad. where 


a = the sum of the lengths of the parallel sides, and 
d 


the perpendicular distance between them. a 


Difficulties arise however, when we consider the area of a region bounded 
by a curve such as a circle or an arch. We can, like Archimedes, cut out 
the appropriate shape and weigh it, and then say that the area is the 
area of the rectangle of the same material and weight. A mathematical 
definition is more complicated and is phrased in terms of limiting 
processes, but it is much more satisfactory. 


In the next section we shall use numerical data to illustrate how we can 
estimate an area (and a volume) in a practical case. The purpose of this 
is twofold. Firstly, we shall see how we can use the areas of rectangles 
to determine upper and lower values (called bounds) between which we 
intuitively know the area to lie. Secondly we shall show that, when actual 
measurements are used in an estimation, there is no need to involve any 
limiting processes, since no increase in accuracy is thus obtained. There- 
fore, unless you are particularly interested in this part of the text, do 
not get too involved with the arithmetic; read through it quite quickly 
and grasp the two ideas mentioned above. 


We shall begin by considering the practical case when part of the boundary 
of the shape whose area we wish to find is specified by a tabulated function, 
where the tabular points are equally spaced. 


We consider a vertical cross-section of a construction site; that is, the 
area we consider represents an imaginary straight slice down through 
the ground. A part of the boundary of sucha cross-section may be specified 
approximately by tabulated values of the heights above a reference 
level of equally spaced points on the boundary. 
We shall see that the accuracy of the estimate of an area which we find 
by taking observations, is usually restricted by the accuracy of the observa- 
tions rather than by the estimation method used. 
We encounter no similar hindrance to accuracy in the calculation of 
area when the bounding curve is the graph of a function such as 
x? 

Six— 10-75 (xe R) 
We shall see that we can then estimate the area as accurately as we wish. 
This will enable us to define the area ina mathematical sense in section 9.2. 


To sum up: in section 9.1 we shall use our knowledge of the areas of 
rectangles to determine other areas as accurately as we can. 


Definition 3 


Discussion 


9.1.1 Region Bounded by a Curve Specified by a Tabulated 91d 
Function 


Discussion 
* 


The sketch above shows a motorway skirting the side of a hill. To con- 
struct this section of the motorway, a large amount of soil and rock 
would have to be removed. To determine the volume of material to be 
transported, any contractor would in practice have to take into account 
the fact that the cross-section of the volume to be removed varies along 
the direction of the road itself, but for illustrative purposes we shall 
consider the problem of determining the volume when the area of the 
cross-section is assumed constant. 


Example | Example I 


As a result of carrying out a survey of the site, the following measure- 
ments were found for the cross-sectional area (of width 20 m) used to 
determine the volume of material to be removed: 


Horizontal distance Height above reference 


in metres level in metres 
0 0 
2 13 
4 3.4 
6 3.9 
8 6.1 
10 7.2 
12 8.5 
14 8.9 
16 10.7 
18 11.8 
20 13.3 


(continued on page 6) 


Solution 1 
Solution 1 


In the above diagram, 
the area of the triangle = 4d(c — b) 
the area of the rectangle = db 
So 
the area of the trapezium = db + 4dc — 4db 
=He+ bd 


= jad | 


(continued from page 5) 


There is a limit to the accuracy with which we can take the measurements. 
Let us assume that there is an error of +0.05m in all the horizontal 
steps and in all the heights. (Cumulative effects of such errors could make 
things worse — this is the best we can do.) The problem is to determine, 
as accurately as we can, the cross-sectional area of the material to be 
removed in this simplified case; the volume can then be determined by 
multiplying by the length of the stretch of road. 


Proposed road 
Surface 


Solution of Example I 


Remember the elementary way we are going to approach this problem. 
We are going to approximate the area by a sum of the areas of rectangles. 
But it would also be desirable to obtain an estimate of the magnitude 
of the error we may be making. We can do this fairly simply by determining 
two areas made up of rectangles: one which will be larger than the area 
we require and one which will be smaller. We then know upper and lower 
bounds of the area we are interested in. We can plot the data on a graph 
like the one given below. 


144 Height above reference 
level in metres 


2 4 6 8 10 12 14 16 18 20 
Put Overlay 1* in position. An extremely crude guess at the answer is to 
say that the area lies somewhere between the area of the rectangle ABCD 
and zero. But we can do better than this. Put Overlay 2 in position. 
You then see an area, made up of rectangles, which is larger than the area 
we require. We have used the heights at 4 m (instead of 2 m) steps in our 
second attempt to estimate the area, in order to demonstrate later the 
improvement we get in accuracy as we increase the number of rectangles 
appropriately. You may well ask at this stage how we know that the 
boundary of the area obtained using rectangles encloses the actual 
surface, since we only have a finite number of observations and there 
might bea sharp peak in the ground in between two of these observations. 
We do not know this, of course, but, as in finite differences, when we have 
the tabulated values ofa function we assume that the graph of the function 
in between the tabulated points is reasonably smooth. (Anyhow, the 


surveyor should have spotted any irregularities.) 


The area of the sum of the rectangles on Overlay 2 is 168.0 m?. We now 
estimate the absolute error bound in this area, due to the absolute error 
bound of 0.05 m in the measurements of height and of the horizontal 
step length. 


We use ¢, to denote the absolute error in the measurement x. In Unit 2, 
Errors and Accuracy, section 2.3, we deduced that 


Cpry = ty + 2, 


and 


Oy = ey + Xl, — B28, 
that is, 

Oxy = Yee Xe, 
since the term ¢,¢, is small. 


* The overlays are in the wallet on the inside of the back cover of this text. 
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Horizontal 
distance 
in metres 


Applying the last result to find the error in the area of each of our rect- 
angles, and then applying the first result to the sum of thesc arcas, we 
get an estimated absolute error in the total area of all the rectangles: 


(0.05(3.4 + 6.1 + 8.5 + 10.7 + 13.3) 
+0054 +44444 +4 4))m? = 3.1m? 


since the absolute error bound in each of our measurements is 0.05 m. 
Thus the estimated absolute error bound in the area is 3.1 m?. 

Now put Overlay 3 on the basic graph. The area shown here is clearly 
less than the area we need (we assume that there are no deep, narrow 
ditches in the ground). 

The sum of the areas of the rectangles is now 114.8 m?, with estimated 
absolute error bound: 


(0.05(3.4 + 6.1 + 8.5 + 10.7) + 0.05(4 + 4 + 4 + 4))m? 
= 2.2m? 


Can we improve this accuracy by using smaller rectangles, and in parti- 
cular, by using rectangles with half the width of the previous ones? 


Put Overlay 4 on the basic graph. You will now see two sets of rectangles 
on a 2m wide base. If you use Overlays 2 and 3 at the same time as 
Overlay 4, you will see that the sum of the areas of the larger rectangles 
with bases 2 m wide is less than the sum of the areas of the larger rectangles 
with bases 4 m wide, but it is still greater than the required area. Similarly 
the sum of the areas of the smaller rectangles with bases 2m wide is 
greater than the corresponding sum for intervals of 4m, but it is still 
less than the required area. 


Carrying out the calculations for this last case we get for the larger area : 
211.3 + 3.44 39+ 61472485489 + 10.7 


+11.8 + 13.3)m? = 150.2 m? 
with an estimated absolute error bound of 4.8 m2* 
Similarly, for the smaller area: 
21.3 + 344+ 394614724854 89 
+10.7 + 11.8)m? = 123.6 m? 


with an estimated absolute error bound of 4.0 m2.* 


Note that, in this example, the difference between the larger and smaller 
areas is simply the area of the largest rectangle. Can you see why this is 
so? This simplifies the calculation, since we can calculate the larger area 
from the smaller area by adding to it the area of the largest rectangle. 
Thus the area has a value between 


(123.6 — 4.0) m? = 119.6 m? 
and : 


(150.2 + 4.8) m? = 155.0 m? 


and we have apparently narrowed the gap between our upper and 
lower estimates by increasing the number of intervals we use on the 
distance-axis from five to ten. The following exercises develop this idea 
a little further. The first of these exercises indicates a method of obtaining 


the “best” estimate ofa number when we know the error interval in which 
it is contained. | 


* Ifyou check these results, remember that we are assumin 


length has an absolute error bound of 0.05 m. oe height and each atep 
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Exercise I 


If you were told that the number of matches in a box was between 44 and 
50 inclusive, how many matches would you say were in a specific box 
if your objective was to make the possible error a minimum? 


What would be the absolute error bound in this case? a 


Exercise 2 


(i) Using the results (obtained on pages 7 and 8) for five intervals in 
our example, and the idea introduced in the last exercise, give an 
estimate of the cross-sectional area in Example 1. State the estimated 
absolute error bound due to the propagation of errors from the inexact 
data, and the absolute error bound arising from the use of rectangular 
approximations. 

(ii) Repeat (i) using the results on page 8 for ten intervals. i 


As we increase the number of rectangles, the estimated absolute error 
bound arising from the inexact data will continue to increase, while 
that arising from the rectangle approximatio. will decrease (as we saw 
in the two parts of the solution to the last exercise). Alter a certain stage, 
this makes the attempt to further improve accuracy by using more and 
more intervals a worthless pursuit. 


Exercise 3 


Describe in words, with the aid of sketches, how you would determine an 
estimate with upper and lower bounds of the volume of material removed 
if the cross-sectional area changed significantly along the line of the road. 


(HINT: Think what to use for your basic element of volume.) n 


Exercise 1 
(2 minutes) 


Exercise 2 
(3 minutes) 


Discussion 


Exercise 3 
(3 minutes) 
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Solution 1 Solution 1 
Forty-seven matches with a possible error of 3. If you choose a value 
other than 47, say 46, your error could be greater, for example, it would be 
4 if there were 50 matches in the box. So the answer is $(50 + 44) matches, 
and the absolute error bound is $(50 — 44). | 
Solution 2 a Solution 2 
(i) Largest possible area = 168.0 + 3.1 = 171.1 m? 

Smallest possible area = 114.8 — 2.2 = 112.6m? 

Estimated area = $171.1 + 112.6) = 141.8 m? 

Estimated absolute 1 2 

= 1 — 112.6) = 29.2 

error bound 171.1 — 112.6) = 29.2m 

of which: 

estimated absolute 

error bound arising = 33.1 + 2.2) = 2.6m? 

from inexact data 

and 

absolute error bound 

arising from the = (168.0 — 114.8) = 26.6 m? 

rectangle approximation 
(ii) Largest possible area = 150.2 + 4.8 = 155.0 m? 


123.6 — 4.0 = 119.6 m? 
(155.0 + 119.6) = 137.3 m? 


Smallest possible area 


Estimated area 


I 


Estimated absolute 
error bound 


of which: 


estimated absolute 
error bound arising = 14.8 + 4.0) = 4.4m? 
from inexact data 


= (155.0 — 119.6) = 17.7 m2 


and 


absolute error bound 
arising from the = 3(150.2 — 123.6) = 13.3 m? 
rectangle approximation 


which is exactly half that in (i), as expected. a 
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Solution 3 ‘ 
Solution 3 


A possible method is outlined below, using a basic element of volume: 


(i) The element of volume is a rectangular prism on a square base. 
We define the volume of such a prism to be the product of the area Definition 1 
of the base and the height. bate 
(ii) A grid of heights is plotted on the surface with equal distances in 
both horizontal directions, for example: 


| 
7.2m 5.3m 3.9m 
61m 5.2m 
1 
(iii) Take the maximum height in each grid square, for example, 6.1 m 
in the one shaded. This height will lead to a rectangular prism based 
on this square with volume greater than the volume required. 

(iv) Take the minimum height in the grid, for example, 3.9 m. This gives 
the smallest volume of a corresponding prism. 

(v) The volume lies between the sums of the minimum and maximum 
volumes of the basic rectangular prisms. The “best” estimate is 
half the sum of the maximum and minimum volumes. 

Alternatively, you could have estimated the cross-sectional areas of 


slices of the site at equal intervals, like slices of a loaf of bread, and then 
found the maximum and minimum values of the volume of each slice. 
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This method has the disadvantage that you need to estimate the magni- 
tudes of all the cross-sectional areas first. 


9.1.2 Region Bounded by a Curve Specified by a Continuous 9.1.2 
Function 


In the previous section we saw that, where there were possible errors in Discussion 
the measurement of a profile, our attempt to continue to improve the ae 
accuracy of our estimated area by increasing the number of rectangles 

used in the calculation was thwarted. It seemed that, if we reduced the 

width of the intervals, we got a better approximation to the area until 

the inherent errors began to dominate. 


{estimate of the area 


the actual 
area f~ 


Number of intervals 


The crosses on the diagram Tepresent a plausible set of calculated values 
for the area (not the actual values) as the number of intervals is increased. 
For a small number of intervals the error is dominated by the error in 
the rectangle approximation, for a large number of intervals by the 
inherent errors in the data. In the example in the last section, the figures 
on page 10 suggest that there would be little point in reducing the 
horizontal steps to much less than | m. In the extreme case, if we can 
only measure the width of the intervals to within 0.05 m, there fg obviously 


not much point in reducing the interval width to less than 0.05 m. (In 
fact, as we have seen, the reduction of the interval width would lose its 
point at a much larger value.) 


If the profile which determines the boundary of the area is given exactly 
at each point, then the last restriction does not apply. This would be the 
case if the profile were the graph of the function 


2 
fix-— 10 — - (x€[—10, 10) 

as in the next example. This is the sort of situation which would arise 
if we were designing a piece of equipment or a building, where we decide 
what shapes we want, as opposed to the Previous example, where the 
shape was provided for us and we had to do the best we could. 

For a region bounded by a curve specified by a continuous function, 
we can, by dividing the domain of the function into as Many intervals as 
we please, approximate the area to whatever accuracy we desire. Does 
this stir memories of the limiting procedures and the concept of continuity 
of a function introduced in Unit 7, Sequences and Limits I? We hope so, 
because this is the way we shall define the area as a limit in section 9.2.1. 
Let us first look at a specific example: 


Example J 


The Dollan Baths, Scotland’s Olympic length swimming baths, were 
officially opened on May 27, 1968, by R. B. McGregor, the Scottish 
International swimmer. 


Designed by architect A. Buchanan Campbell, the building rises in 
one immense parabolic arch, which has a span of 100 metres and which 
reaches a maximum height of 20 metres above ground level. 

When he had decided on the shape, the architect had to calculate the 
cross-sectional area of the building to enable him to find the pressure 
exerted on the structure. 

In order to emulate his calculation, we first need to find a function f 
whose graph is a curve which represents the outline of the building on 
a suitable scale. 
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Example 1 


You may like to try to find this function, knowing that its graph is part 
ofa parabola, which means that it is a function of the form 


x-—a+bx+cx? — (a,b,ce R) 
| 


x fame 
You should get the answer x-—> 20 — De However, our aim is to 


calculate the area underneath this graph. There are two parts to this: 


(i) we approximate to the area with rectangles, find upper: and lower 
estimates of the cross-sectional area of the parabola and hence find 
the best estimate of the area in the two cases: 

(a) when OC is divided into five equal intervals; 
(b) when OC is divided into ten equal intervals; 

(ii) we determine how many intervals of OC would be required to ensure 

that the estimated absolute error bound is less than 1 m?, assuming 


that exact arithmetic is used. 


Solution of Example | 


(i) Have a look at the basic graph without any of the overlays* in place. 
The parabolic arch is symmetric; therefore, the area represented by 
AOCB is twice the area represented by OBC. The images for the given 
values of the variable x are calculated exactly: 


0 | 5 30 35 | 40 45 | 50 
a ae 15 | 12.8 102 | [72 38 | 0 


(a) Now put on Overlay 5. 


x 
f(x) 


Using five intervals for the half-area, form the area consisting of 
rectangles which give an area larger than the area required. 
These rectangles are outlined by solid lines in the figure. This 
gives us a first upper estimate for the total area: 


2 x 10{20 + 19.2 + 16.8 + 12.8 + 7.2} m? = 1520 m? 


The dashed lines in the figure indicate an area made up of rectangles 
which give an area smaller than the area we require. 


This gives a first lower estimate for the total area: 
2 x 10{19.2 + 168 + 1284724 0} m? = 1120 m? 


We notice that the difference between the two estimates is twice 
the area of the largest rectangle. 


* The overlays are in the wallet on the inside of the back cover of this text, 


Thus the best estimate for the total area using five intervals is 


1520 + 1120 


5) m? = 1320 m? 
with estimated absolute error bound 


1520 — 1120 


5 m? = 200 m? 


(b) Remove Overlay 5 and put on Overlay 6. 


When ten intervals are used for the half area, we find that the 
upper bound on the area using ten rectangles is now 


2 x 5{20 + 19.8 +++. 47.2 + 3.8} m? = 1430 m? 


Calculation of the lower bound gives an area = 1230 m? 


Our estimate would be pa m? = 1330m? 
with estimated absolute error bound now = 100m? 


(ii) We are now looking at a different aspect of the problem. 


We state at the beginning how accurate we want our estimate of the 
total area to be. We want it to be accurate to within +1 m?, that is, 
we require an absolute error bound of 1 m?. 


For what width of interval can we achieve this? As we have already 
noticed, the difference between the upper and lower bounds of the 
estimate of the area is twice the area of the largest rectangle, and it 
is also twice the estimated absolute error bound. Therefore the area 
of this rectangle must be less than or equal to 1 m?. But in this example 
we know that the height of the largest rectangle is equal to 20m 
regardless of the number of intervals. Therefore, to achieve the 
desired accuracy, its width must be at most 0.05 m, since 


20 x 0.05 m? = 1 m? 


The total number of intervals along OC = aa = 1000. So 
interval width 


to be sure of achieving the desired accuracy by this method, we must 
divide OC into at least 1000 intervals. 

Similarly, we can show that, to achieve an estimate of the total area 
with absolute error bound 0.1 m?, we would require at least 10 000 
intervals; and that in general to achieve an estimate with absolute 
error bound em? we would require at least n sub-intervals where n 
is an integer such that 


1000 
g 


n2 


\ 


N.B. Later in this text we shall see that we can find the above parabolic 
area exactly. Its value is 1333}m7. So, in fact, the above estimates are 
very good ones with an error which is already less than 15 m? in the 
five-interval case. | 


Exercise 1 


In the two examples discussed in the text, we remarked that the difference 
between our upper and lower estimates of the required area (using 
rectangles) was simply the area of the largest rectangle (see pages 8 and 
14), Can we say this when trying to find the area under any curve by 
using rectangles, or was there something special about these two 
examples? i | 


Exercise 1 
(5 minutes) 


Solution J 


The statement we made applies in the general case only if the graph of 
the function slopes downwards or upwards over the whole region. 
example, the statement is not true for the diagram shown below, althoug! 
we can still use rectangle approximations to find the area. One set of 
rectangles could be chosen, the sum of whose areas is greater than the 
required area, and the other set so that the sum of the areas of these 
rectangles is smaller than the required area. An example of the first set 
of diagrams is shown in the figure. 


Ric) 


Xref (x) 
x 
| 
Exercise 2 
treo 
20m4 
nvm 200 

10m-4 

> 

10m 50m * 


The boundary of the cross-section of the roof of a main airport building 
is the graph of the function 


fie (x € [10, 50]) 


The supporting walls are vertical and situated at x = 10 mand x = 50m 
respectively. 


(i) Using four intervals, estimate the cross-sectional area (indicated in 
the diagram) as accurately as you can by the method we have used 
in the example, and state the accuracy which you think you have 
attained. 


(ii) How many intervals would you require to attain an estimate with 
absolute error bound 1 m2? a 


Solution I 


Exercise 2 
(5 minutes) 


Exercise 3 Exercise 3 


(4 minutes) 


X+>COS X 


Nia 
~ 


Given the function 


fix'— cos x (<[03]} 


we can find an estimate for the area between its graph and the two axes, 
by forming rectangles in the following way: 

We divide the domain into equal intervals and take the height of the 
rectangle on each interval to be 4 (the sum of the ordinates* of the end 
points of the interval). E.g. the height of the rectangle on EA = (AB + EC). 


Are the following statements true or false? 


(i) The sum of the areas of the rectangles above is 

the same as the estimate of the area which we have 

used so far. That is, the estimated area is the same 

as that obtained by the previous method: } (the 

sum of the areas of the larger rectangles + the sum 

of the areas of the smaller rectangles). TRUE/FALSE 
(ii) The previous method of estimating the area is 

better than the method described in this exercise 

because it gives us an estimate of the magnitude 

of the error. TRUE/FALSE 


* The ordinate of a point in the Cartesian plane is the y-co-ordinate of the point. 


Solution 2 


(i) x 10 | 20 30 | 40 | 50 


4 


f(x) | 20 | 10 667 | 5 


The total area of the larger rectangles = 416.7 m2. 


(Notice that in this exercise an inherent error, due to the finite decimal 
representation of 3, has been introduced. However, we can always 
make this as small as we please.) 


The total area of the smaller rectangles = 256.7 m?. The estimate for 

416.7 + 256.7 . 
the required area is therefore (67 + 2567) 2 = 336.7 m?, with 
estimated absolute error = a m? = 80 m?, so in this case the relative 
error is large. 


(ii) The differences in height and area between the largest and smallest 
rectangles are 16 m and 16 x (interval width) m? respectively, so we 
require 16 x (interval width) = 1. Therefore the width of each 
rectangle must be gm. The total number of intervals required is 


therefore a = 640. 
16 


We noticed previously in Example 1 that the difference between the sum 
of the larger rectangles and the sum of the smaller rectangles is a very 
crude estimate of the error. In fact, the exact area in this case can be 
calculated and it is 321.9 m? (to one place of decimals), which lies well 
inside the error interval which we obtained in (i). | 


Solution 3 


(i) TRUE. When in our original method we take the estimate } (the sum 
of the areas of the larger rectangles + the sum of the areas of the smaller 
Tectangles) this is exactly the same as the sum of $ (the area of the 
larger rectangle + the area of the smaller rectangle) taken over all 
the intervals. 

(ii) TRUE. The original method also told us the magnitude of the error 
we might be making. Although this error estimate may be a wild 
exaggeration, it is better to have some estimate than none at all. 


Solution 2 


Solution 3 


9.1.3 Summary 


In this section we have explored the idea of estimating the (intuitive) 
area of a region bounded by a curve. The methods we have discussed 
have considerable practical application, to which we return in section 9.4. 
For the moment, however, we wish to use the ideas of this section to 
develop a mathematical definition of area under a curve. 


9.2 THE DEFINITE INTEGRAL 
9.2.0 Introduction 


In section 9.1 we developed a method of approximating to what we 
intuitively call “area”, by considering a number of rectangles, this number 
being determined by the accuracy required, except in the cases when 
inherent errors occur. In this section, we shall tie up the intuitive definition 
of area with a mathematical one based on the ideas you met in section 9.1. 
We assume that you were satisfied with the statement that, in the absence 
of inherent errors, increasing the number of rectangles increases the 
accuracy of the estimate of the area. If so, you were quite justified, but 
in order to help you to see that we were making an assumption, let us 
consider the problem of estimating the length of a curve. 


Exercise | 


We consider a line AB of known length to find out if an estimation pro- 
cedure, similar to the one we have been considering for areas, gives us a 
reasonable approximation to the known length. 


The graph of the function 
x-—1-x (x € (0, 1)) 


is a straight line, AB. 

We approximate the length of this straight line by a “stair-case”’, the 

stair treads (intervals) being of equal length. Just as in the area estimation, 

we can do this in a number of ways, of which we illustrate one. 

(i) In the sense of finding the length of a stair carpet, find the total 

length of the zig-zag line in the figure above. 

(ii) Does the number of steps in a stair-case like the one illustrated make 
any difference to its total length? 

(iii) Intuitively does it appear to you that we can get the corners of the 
stair-case as near to the line AB as we like by taking enough steps? 
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Summary 


9.2 


9.2.0 


Introduction 


Exercise 1 
(2 minutes) 


(iv) Do you think that this implies that we can take enough steps to 
make the length of the staircase approximate as closely to the length 
of the straight line as we want? 

(v) What then is the length of the line, in the limit, by this procedure? 

(vi) Is this the correct length of the line? If not, what is wrong? a 


The purpose of this last exercise was not particularly to investigate the 
idea of length but to show that we must be careful in the limiting pro- 
cedures we adopt to back up intuitive ideas. Approximating a curve by a 
“staircase” seems satisfactory when finding the area bounded by the 
curve, but it is unsatisfactory when trying to find the length of the curve. 
In both cases we are trying to generalize a concept. Whenever we do this, 
we must check that our generalized definition gives the answer we expect 
in the more fundamental case. The ‘‘staircase” approach does not give 
us the “right” answer for the length ofa straight line: so there is no point 
in trying to generalize it. But our method for the calculation of area 
does satisfy intuition and accords with our concept of area in simple 
cases. Further, we were able to sandwich the required value between two 
estimated upper and lower bounds. Thus we can speak of the “‘accuracy” 
of our estimate, since we have trapped the numerical value of the area 
in an error interval which we can make as small as we please, and so we 
can proceed with confidence. 


In section 9.2.1 we shall extract the mathematical definition of the definite 
integral which we have chosen to use. We say “‘chosen to use” because 
there are several ways of defining the definite integral, some of which 
are more general than others. 


In section 9.2.2 we shall show that it is not always necessary to determine 
definite integrals completely from scratch. 


We find that functions such as: 


xe— 3x — x? +2 (xe R) 
x 2sin x — 2 (xe R*) 


and other similar functions can be expressed in terms of simpler functions, 
such as: 


x—> 1, xe x, x x?, x x3, x sin x 
(xe R), 
1 1 
oS (xe R, x # 0) 


(See Part 2 of Unit 7, Functions, in which we discussed combinations of 
functions.) Can we find standard expressions for the definite integrals of 
simple functions such as these? If so, can we use them to find the definite 
integrals of more complicated functions, and thus simplify our work? 
The answer to the first question is “yes”. We proceed to find some of these 
answers in section 9.2.1: in other cases the calculations would be too 
lengthy at this stage, and so we shall leave them to Unit 13, Integration II 
The answer to the second question is also “yes” in many cases, such as 
the examples given above. For combinations of simple functions such as 


xr—>3x—-—x?4+2 (xe R) 


we need to develop the rules for combination of the appropriate definite 
integrals, and this we do in section 9.2.2. 
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Discussion 


9.2.1 Area and the Definite Integral 


Using the experience we have gained in section 9.1, let us try to extract 


the definition of a definite integral in such a way that it matches up with 
our intuitive idea of area. 


In section 9.1 we took a region and sandwiched it between two sets of 
rectangles, one containing the region and the other contained in the 
region. We agreed that the “‘area” of the region lay between two areas: 
the sums of the areas of each of the two sets of rectangles. 


Suppose that we use this method to obtain upper and lower estimates 
of the area of the region by considering two sets of rectangles, each with 
n members. 

We denote the sum of the areas of the larger rectangles by ,,and the sum 
of the areas of the smaller rectangles by «,, .1 can be any positive integer, 


so if we let n take the values 1,2,3,... successively, we obtain two 
sequences: 


and 
Oypetlyo Base es 


which we denote by 1 and @ respectively. (See the notation for infinite 
sequences introduced in Unit 7.) 


The terms of A are always greater than the required area, and the terms 
of g are always less than the required area, but as n increases, the corre- 
sponding terms in the two sequences (intuitively) get closer and closer. 
So as the number of rectangles gets very large, we can say (intuitively) 
that the two sequences have the same limit, and we can define this limit 
to be the arca of the region. We have to be a bit more specific about area 
when we think in terms of more general functions. Let f be a function 
with domain [a, b], such that the graph of f, the lines specified by x = a 
and x = b, and the x-axis form the boundary of a closed region. 


‘Area beneath the 
graph ot t, 
between a and b. 


The area of this region is usually called the arca beneath the graph of f 
between uw and b even though the graph of f may look like this: 
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9.2.1 


Definitions and Notation 


Notation 1 


Notation 2 


Definition 1 


Definition 2 


(continued on page 22) 


Solution 9.2.0.1 Solution 9.2.0.1 


(i) Each rising tread is equal to a corresponding length on the (1 — x)- 
axis. These lengths do not overlap, and therefore they add up to the 
length OB. 

Each step length is equal to a corresponding length on the x-axis, 
and again the lengths do not overlap, but add up to the length OA. 
So 


the total length of the zig-zag line = the total length of the 
carpet required 


=1l+1=2 


(ii) NO. The total length remains unchanged no matter how many steps 
are used to cover the line. 

(iii) Intuitively YES. 

(iv) You have probably answered YES. 

(v) Apparently 2 units, using (i) and (ii). 

(vi) NO. The correct length is J2 units. The stair-carpet, or zig-zag line, 
is always either horizontal or vertical and is never lying flat on the 
slope of the hypotenuse of the triangle. So, although, the answer to 
(iii) is correct, it is misleading. | | 


Se 
(continued from page 21) 


or this: 


At the end of the example about the Dollan swimming baths (page 15) 
we found that to achieve an estimate of area with absolute error bound 
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em? we would require at least n intervals, where n is an integer and 


Ace 1000 
é 


\ 


That is, we found we could choose n such that the estimate of the area 


belonged to an error interval of width at most 2e, where é could be as 
small as we pleased. Here we have a similar case. 


Expressed the other way round, this means that as n increases the width 
of the error interval decreases. 


An 


Actual 
Area . eG 


an 


For n intervals, the upper and lower bounds of the error interval are 

determined by A, and a, respectively. 

If lim A exists, then lim g will also exist, and vice versa. The two limits 

will be the same number: call it A. As n increases, the width* of the error Notation 3 
interval will shrink towards zero, and the area of the parabolic section, 

which is.always contained in the error interval, must (intuitively) therefore 

also be A. 


Now let us consider the problem for the function f more algebraically, 
and pick a particular sum (you will see why in a moment): 


S, = hf(a) + hf(a +h) +--+ + Af(a + {n — 1h) Definition 3 
= hl f(a) + fla +h) +--+ + f(a + {n — 1}h)] 


where hh = tat that is, we have divided [a, b] into n sub-intervalst each 
n 


of width h. This expression for S,, defines a sequence S,.S,,53,... which 
we denote by S. For a particular function f with a graph as illustrated : Notation 4 


* The width of [a, 6) is b — a. 
+A sub-interval of [a, b] is an interval [c,d] wherea < ¢ <d <b. 


we would have 
S, = Any 


the sum of the areas of the larger rectangles. If f has a graph like this: 


h= a 
then 
S, = an 
If f has a graph like this: 
ies) 
Xre f(x) 
Hi 
a b x 
then 
a, < S, <A, 


The point is this: if the graph of f isa smooth curve, and f(x) is positive* 
for all values of x in [a, b], then, provided that 


limg=lim4=A 
it follows that 
lim $ = A 


since S, will be sandwiched in an error interval which diminishes in size 
as n gets larger (just like the estimate of the parabolic area in the Dollan 
Baths example). 

Because A, the limit of S, has tremendously important uses in contexts 
other than that of area, we give it a special name. We call it the detinite 


integral of the function f between a and h (or in {u. bj), and a special 
symbol 

" 

fi 


Ju 


* The reason for requiring the images to be 


' Positive at this stage will become clear from a 
consideration of Exercise 1 which follows. 
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Definition 4 


Notation 5 


so that 


The symbol [ an elongated ‘‘s”, represents the summation process with 


the “a”? showing where the sum begins and the ‘b” showing where it 
finishes. (The “a” and “‘b” are called the end-points* of the integral). 

We can now forget about the rectangles, and concentrate on this par- 
ticular sequence of sums, S, and its limit. We do this because we are 
interested not only in area but also in the importance of the definite 
integral in mathematics and applications. Sometimes the integral will 
represent an area, but on other occasions it may represent a volume, or 
an average, or an electric current, or a probability, or power, or distance, 
or a similar quantity. In each case, the definite integral will have to be 
interpreted with caution. To avoid some pitfalls we shall stipulate one 


rb 
condition on f which guarantees the existence of f Jf. Remember that 
a 


not all sequences converge (i.e. have a limit) and we have defined the 
definite integral to be the limit of a sequence. We do not want to have to 
return to the definition repeatedly to check on the existence of integrals. 


b 
So we ask: Are there any special types of function f for which f S always 


exists? The answer is “Yes”, and in fact we are already on solid ground. In 
our definition on page 24 we have included the phrase “‘if the graph of fis 
a smooth curve”’; that is, the graph of f must have no “gaps” in [a, b]. 
This is the same as saying that “the function f must be continuous in 
(a, b}” (see Unit 7, Sequences and Limits I). It can be shown (by a more 
rigorous treatment than ours) that if f is continuous in [a, b], then the 


b 
definite integral f automatically exists (i.e. its existence is guaranteed). 


a 
We ask you to accept this fact. (It means that the words “provided that” 
before “lim @ = lim A = A” on page 24 are redundant and may go 
home.) Hereafter, in this text, we shall assume that the definite integral 
exists (i.e. the sequence $ converges) in all our discussions. 


Summary 


Given an interval [a,b] which is divided into n equal sub-intervals of 
width 

_b-a 

OR 


h 


and a function f which is continuous in [a,b], we define the definite 
integral of f in [a, b] as the limit of the sequence §, where 


n= At f(a) + fla +h) +--+ + fla + tn — Ih] 


We write 


lim $ = [s 


Exercise | 

Is “the definite integral of f between a and b” synonymous with “the 
area beneath the graph of f between a and b”? 

HINT: Consider the two definitions for the diagrams on page 21. | | 


* Some authors use the word limtis. 


Definition 5 


Summary 


Exercise 1 
(3 minutes) 


Solution | 


Your solution should contain the following points: 

(i) Even if the curve is wholly above the x-axis, ifit goes up and down then 
the two sequences A and $ are not the same, but, intuitively, the terms 
of the sequence § lie between the terms of the two sequences @ and A. 
If these latter two sequences have the same limit, then so also has § 
and we may therefore give an intuitive YES. 

(ii) If some of the curve lies below the x-axis, then we have negative 
terms in §, hence lim $ is definitely not the area. So the definite integral 
and the area may both exist and be different. | 


We may, therefore, conclude that as long as we are careful to ascertain 
exactly where the curve representing f lies, we can use the definite integral 
to find the area under the graph of f between a and b. If f(x) is not positive 
for all x € [c, d], where [c, d]is a sub-interval of [a. b] (see diagram), then we 
find the definite integral of f between 


fac], [ed], — [d,b] 
separately, and adjust the sign of the definite integral of f between c and 
d before adding the three results to find the total area. 


‘ 
f(x) 


The following exercise is designed to give you practice in the ideas intro- 
duced in this section and to check that they are in accordance with 
intuition. 


Exercise 2 
Find the definite integral of the function f in (a, b], where 
fixe 1 (xe [a, b], b > a > 0) 


and check whether it does give the area beneath the graph of f. nH 


Exercise 3 
Find the definite integral of the function J in [a, 6), where 
fixt— x (xe [a,b], b > a> 0) 


and check whether i i 
be ‘a whether it does give the area beneath the graph of f between 


HINT: In the solution you will need the fact demonstrated in Unit 4. 
Finite Differences, page 47, that the sum of the first n natural numbers is 


Sis Lee 1) . 


Solution 1 


Discussion 


Exercise 2 
(S minutes) 


Exercise 3 
(3 minutes) 


Just a word on notation. We have written the definite integral ofa function 
f in [a, b] as 


b 
fr 
a 
In particular cases when f is known, for example: 


fix x (x € [a, b]) 


we write 
b b 
i) x-— x, or I (x--— x) 
a a 


and omit the domain of f because the part of the domain we are interested 
in is clear from the “a” and “b” in the notation. 


In many text books you will find the notation 


[10 dx 


or, in particular, 


b 
f xdx 
a 


This is a notation which has been used for many years. If youare interested 
in its origins, you will find it in a book on the history of the subject. The 
reasons for our choice of notation are discussed in the introduction at 
the beginning of this text. 

The definition of definite integral we have used is not the most general 
one. We have, for example, chosen the particular case where we divide 
{a, b] into equal intervals. The definite integral of f in [a, b] can be satis- 
factorily defined in a more general way, omitting the conditions: “‘f is 
continuous in (a, b]” and “the sub-intervals of [a,b] are of equal width”. 
On the whole we have relied more on intuition rather than mathematical 
rigour ; we leave the rigorous definition of the definite integral to a later 
course. 


In the last two exercises you found that 


bo, b b 
foren-b-a and fo en=5* 
b>a>0 


By a similar method to that used in the exercises, it can be shown that 


b Ba 
(xr x?) = a 


a 


To derive this result, one needs the formula for the sum of the squares 
of the first n natural numbers 


= n(n + 1)(2n + 1) 


S2(n) 6 


which can be found in Unit 4, Finite Differences, page 48. 


Notation 1 


Main Text 


(continued on page 29) 


Solution 2 


The interval is divided into n sub-intervals of width h, where nh = b — a. 
Then 
S, = hALf(a) + fla +h) +--+ + f(a t+ {n - 1h)] 
=h[1+1+---+1)=nh=b-a 
(n terms) 

The definite integral between a and b is the limit of the sequence §, that 
is (b — a). 

This is clearly the area beneath the graph between a and b. a 


Solution 3 


The interval is divided into n sub-intervals of width h, where nh = b — a. 
We have f(x) = x, so that 


f(a) =a,f(a+h)=a+h,...,f(a + {n- Ih) =at {n— Ih 
Therefore 
S, = AL f(a) + f(a + h) + f(a + 2h) +--+ f(a + {n — A)] 
=h{a + (a +h) + (at 2h) +--+ + (a+ {n—- Ih} 


=nha + WIL + 2434-4 {n- 1] 
(collecting together like terms) 


_— 7)? 
ab — a) +2 — 4) (L+2+34---4+{n-1}J 


(substituting for h) 


i] 


a ie 
alb —ay+2 9 Sin — 1) 


(b — a? (n-1)n 


bess 
a a) + 2 5 


il 


(substituting for $,(n — 1)) 


Solution 2 


Solution 3 


The limit of this sequence for n large is 


b — a2 Cee | 
a(b a+! me a = 2) [sinc lim (:-4)=1] 


2 


nlarge 
So the definite integral of f in (a, b) is 4{b? — a”). 
By regarding the area beneath the graph of f between a and b as the 
difference between the areas of the larger triangle, whose base is the 
X-axis from the origin to’b, and the smaller triangle, whose base is the 
x-axis from the origin to a, we see that this area is 
4? — da? 


which is the same as the definite integral of f in [a, b]. = 


(continued from page 27) 


Exercise 4 Exercise 4 


. . F - (2 minutes) 
(i) From the results so far obtained, an intuitive guess at the value of 


b 
f (x-— x"), b>a>0, meZt 
a 


might be 


pnt m+ 


(a) — 


m 


bm — gmt 


m 


(b) 


prt _ gmet 


©) m+ 


bn — g™ 
m 


(d) 


Which do you think is correct? 
(ii) To obtain the above value would you need to know the sum of the 


(a) (m — 1)th 

(b) mth 

(c) (m+ 1)th 

powers of the first n natural numbers? | 


Exercise 5 Exercise 5 
(2 minutes) 
Calculate, using the general result for the definite integral of x-—— x”, 


obtained as the correct answer to the previous exercise, the values of 


4 

(i) J (x-— x3) 
2 

(ii) J (xr x4) 
-2 


1 
(iii) ik (x-— x43) a 


Solution 4 Solution 4 


(i) (c) 

(ii) (b) | 

Solution 5 Solution 5 
(i) 64 

ji) 

(iii) a5 a 

You have probably recognized that the calculation of Discussion 


b 
f (x-— x") (me Z*) 


from first principles becomes more cumbersome as m increases in magni- 
tude. That is why we leave it until we get some tools in Unit /3, Integration 
II, which can handle it very simply. The formula we obtain if we did the 
calculation is, as already noted, 


mt) amt 


a 


+ 
m+i1 (mezZ*) 


b 
I (x-— x") = 
a 
In fact, although we cannot show it at this stage, this formula is true for 


(me R, m # —1), and you may assume this for any subsequent exercises 
in this unit. 


Table of Definite Integrals of 
Simple Polynomial Functions 


b 


a 
xr— I] b-a 
<7 Xx b? = a? 
xt x? bP =a? 
3 
pet gmt 


Xt— x™ (me R, m# —1) ra 
m 


The first two above were derived in the text; we indicated how to prove 


the third, and you may use the general result although we have not proved 
it. 
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9.2.2 The Definite Integral for Combinations of Functions 
In order to extend the class of functions for which we can calculate 
definite integrals to include, for example, 
| xt-—+ (3x3 + 5x) 
1 


we need two theorems (unless, of course, we resort to summing series!). 
We derive one of them in the following example and the other in the 
exercise below. 


Example 1 


Show that 


[ Wee ef's 


where f is a function, c is any number and ¢f is the function 


xr of (x) nm 


Solution of Example 1 


This is not sucha rigorous proofas that which appears in B. Hunt, Calculus 
and Linear Algebra (W. H. Freeman 1967), page 273 (see Bibliography), 
but the essence of the argument is here. 


To define the definite integral of f between a and b we used the limit of 
the sequence §, where 


S, = hE f(a) + fla +h) +--+ + fla + {n — 1h) 


Similarly, to define the definite integral of (cf) between a and b we would 
use the limit of the sequence T, where 


T, = h{cf(a) + of(a +h) +--+ + of(a + {n — Ih) 
= cS, 


It follows from the result on page 33 of Unit 7, Sequences and Limit 1, 
that 


lim T = clim $ 


and hence 
b b 
[anefs " 


Exercise | 


Convince yourself, by using the definition of the definite integral, that 
if f(x) and g(x) are positive for all values of x in [a, b], then 


[ (f +8) = [s + [fs a 


In section 7.3.1 of Unit 7, Sequences and Limits I, page 34, we intuitively 
derived the result 


lim (§ + T) = lim § + lim T 


(A more rigorous proof was given in an appendix.) In that section of 
Unit 7 we interpreted this result in terms of morphisms, and we may 
similarly interpret our present result. 


al 


9.2.2 


Main Text 


Example 


Exercise 1 
(3 minutes) 


(continued on page 32) 


Solution | 


Write 


S, = hl f(a) + flat hy) +--+ + fla + fn — 1h) 
T, = lela) +g (ath) +--+ + ea + fn — fh] 


and take the appropriate limits. a 


(continued from page 31) 
Let F be the set of all functions f for which it is possible to evaluate 
po 
[7s 
Ifa and b are fixed real numbers, then the function J, defined by 
b 
rp—f se fer 
a 


maps the elements of F to the set of real numbers. 
The result of this exercise then becomes 


Kf + g) = I(f) + Ug) 


where the + on the left represents the addition of functions, and the 
+ on the right the addition of real numbers. We can express this result 
by saying that J is a morphism of (F, +) to (R, +). 

The appropriate commutative diagram is 


(f,8) —*——> ft+eg 
1 i| 
(US), Hg) ——— If) + Ig) = IU + g) 


Noticing that, for instance, 


1 
[x—ciet 


° 


and 


1 
i) Xh— 2x = 1, 
0 


we see that I is many-one, and hence it is a homomorphism. In fact, 
we could have approached this section from this point of view, first 
noticing the function J and then investigating whether it is compatible 
with any of the known binary operations on F, for example, +. x and o 


(the composition of functions). The answers for x and e is that it is not 
compatible. 
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Solution 1 
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Example 2 
Example 2 
Calculate 
a 
xXt— (3x3 + 5x 
f ( x) = 
Solution of Example 2 
3 
(x-— (3x3 + 5x)) 
1 
3 3 
= f (x-— 3x3) + ) (x!— 5x), from rule in Exercise 1 
1 
3 3 
=3 ih (x x9) + sf (x+— x), from rule in Example 1 
1 
34 _ 14 32 _ 42 
=3 
cpus 
= 80 | 
Exercise 2 Exercise 2 
Calculate Giminutes) 
4 
(i) i x—> (2x? + 7x — 3) 
2 
2 
Gf waxy 
-2 
3 
(iii) i x—> (x? — x) 
2 
3 
(iv) i XH—>(x — 1)(x — 2) 
0 
(v) What is the area beneath the graph of f in (iv)? a 
It is frequently possible to find areas of more complex regions with very Main Text 
little extra effort. ah 
Suppose we wish to find areas which are bounded by curves and lines 
specified by equations of the form x = c, ce R, but are not bounded by 
the x-axis ; for example, the area shaded in red in the diagram. 
#(x),g(x) aaitsl 
x2 (x) 
x (continued on page 35) 
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Solution 2 


4 
{i) i? (x (2x? + 7x — 3)) 


= 2f ox) + 7f 9 - af eo 


3 2_ 52 
=o" 57) (8 2 34 — 2) 


2 


3 


_ 220 
3 


The graph of the function is shown in the diagram, and because part 
of the graph of this function lies below the x-axis, the definite integral 
evaluated in (iv) does not give the required area. So we divide the definite 
integral into 3 separate limiting sequences. In the first and third sequence 
all the terms are positive whilst in the second they are negative. The 
appropriate definite integrals are 


1 
i (x (x? — 3x + 2) = 2 
0 


2 
f (xt (x? — 3x + 2)) = -$ 
1 


3 
I xr (x? = 3x + 2) = 8 
2 


The area is the sum of the magnitudes of these three definite integrals 
a8 ed 8 
=ErerEr Ss: a 
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Solution 2 


Rather than try to form rectangles to determine the numerical value for 
this area, we can consider it as the difference of two areas bounded by 


the x-axis. Look more closely at part (iii) of the last exercise from a 
graphical point of view. 


3a 


Since the images of x under the functions x-—> (x? — x), xe—— x7, 
x'— x are not negative for x in (2, 3], the definite integral 


f x-—> (x? — x) 


represents the area shaded in red in diagram (a), and the definite integrals 


f (x x?), i (x-— x) 


= 2 


represent the areas ABEF and ABCD respectively in diagram (b). 
The shaded areas in (a) and (b) may not look equal, but since 


[io ot =a) = foe x) fie > x) 


these areas are the same. 


ts 
a 


(continued from page 33) 


Exercise 3 


Write down the definite integral, or sum of definite integrals, which 
represents the areas on the following graphs, and hence evaluate them. 


(i) sfx) 


A 


x8 


(ii) 


(iii) 


For this one, you will need the result 


b 
f x > x2 — 2(p3/2 “93/2 
a 


(iv) Is it a coincidence that (ii) and (iii) have the same answer? | 
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Exercise 3 
(5 minutes) 


Exercise 4 


Exercise 4 
4 (3 minutes) 


f(x) 


If the shaded area represents a piece of metal 1 cm thick, and x and its 
images are measured in centimetres, calculate the volume of material in 
such a metal plate. a 
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Solution 3 


1 2 
(ii) x4 (8x — 2°) + f xt-—> (8 — 2x?) = 22 
oO ‘3 


The first definite integral represents the area shaded in red, the 
second the area shaded in black. First find the appropriate inter- 
sections of the graphs, then write down the integrals to give the result 
above. 


. I 
Wa xe— |— x x2 — 
(iii) I BR 
(iv) NO. In (iii) we have drawn the graphs of the inverses of the functions 
used in (ii). Effectively we could have used the same diagram as 
in (ii) and interchanged the axes. We often face the choice of forming 
the definite integral “the other way”, that is, using the inverse 
functions. The simplicity gained in the representation of the area 
can, however, be lost in the awkwardness of the new function, as 
we see in this example. | 


fe 


Solution 4 


38 


Solution 3 


Solution 4 


If the areas above and below the x-axis are calculated separately by 
definite integrals, we have 


2 2 
f (xt x?) + the magnitude of [I (xe -| 
0 0 
=$+2=4 


so the volume is 4cm>. We can also represent the whole area by the one 
definite integral 


2 ° 4 
[a —+ (x? — (—x)) [» +(x? + x) = 4, 


fey °O 
and again we find that the volume is 4 cm3. a 


9.3. SOME APPLICATIONS OF THE DEFINITE 
INTEGRAL 


9.3.0 Introduction 


We have defined the definite integral of f in [a,b] as the limit of the 
sequence 


S1,S2,S3,-005 Shy: 
where 
S, = hf (a) + fla +h) +--+ fla + {n — Ih) 


We have also shown that provided the images are positive, this integral 
will give us the area beneath the graph of f between a and b. In this section 
we intend to look at other physical situations which, when we analyze 
them, produce sums of terms which we can identify with a definite integral. 
For example, the terms in the sequence S,,S2,53,...,5,,-.., which 
were used to represent the sums of the areas of rectangles, can be used 
to represent the sums of volumes of standard shapes. We look at this in 
more detail in the next section. 


39 


9.3 


9.3.0 


Introduction 


9.3.1 Volume of a Solid of Revolution 


Consider a region bounded by the graph of a function f (whose images 
are positive in [a,b]), the lines specified by x = a and x = b, and the 
interval of the x-axis between a and b. Focus attention on the boundary 
of such a region, and imagine the boundary being rotated round the 
x-axis; so that it generates the bounding surface of a solid. We usually 
refer to the rotation of the graph of f in [a, b] (leaving the rotation of the 
area to be understood.) The volume of such a solid is called a Volume of 
revolution, and the calculation of these volumes is simply an application 
of the definite integral. 

As an illustration, we shall investigate the problem of finding the volume 
of a cone generated by rotating the hypotenuse of a right-angled triangle 
around the x-axis. 


You may know that the volume, V, is given by the formula V = $nR?H 
where H is the height of the cone and R is the radius of the base. We shall 
derive this formula. 


To correspond to our elementary rectangles (used in the calculation of 


area) we choose discs of thickness h. They are similar to the discs on the 
child’s toy shown in the diagram. 


If there are n discs we have 
H 


h=—, 


n 
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9.3.1 


Main Text 


Definition 1 


and the edge of the cone will be the graph of the function 


Rx 
x (xe (0, HI) 


Rx 
H 


The radius of each of the n discs will be given by an ordinate. Let the n 
ordinates be at x9,x,....,X,—}, then 


Xm = mh (m= 0,1,2,...,2— 1) 


and the corresponding ordinate is the image of x,,, so that 


and the volume of the elementary disc is therefore 
J Rmh\? h 
H 


The total volume of these n discs is 


I 
V,=0+—; P+ am x mee. 
R2V3 
+o x (nF 
R? R? R? 
ifo + IP be + phy? +o + fin - op | 


We write the expression for V, this way because we are going to compare 
it with the sum whose limit we know to be a definite integral, rather than 
go through the lengthy algebraic process of finding the sum again. Thus, 
comparing V, with S,, given on page 23: 


S, = ALf(a) + f(a +h) +++» + f(a + mh) +--+ 
+f(a + {n — 1}h)} 


it is fairly easy to see that the appropriate function is 
-R? 
fixe sara 
and that the end-points of integration area = Oandb = nh +a=nh=H. 
Thus 
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Old hands at the game would abbreviate the above comparison even 
further. They would jump straight from the volume of the elementary 
disc in the mth position: 


Rx,,\? ; 
estat 
a : 
to the appropriate function and then to the integral, by dropping the 


“h” and the subscript ‘tm’. In fact, the justification for doing this is 
outlined by the fuller argument above. 


We can generalize the definite integral for finding volumes of revolution. 
If the volume of revolution is formed by the graph of f rotating about the 
x-axis between x = a and x = b, then the volume generated is 


b b 
f —eatsony? or xf x {sta} 


b 
(We could also use nf f? but f?, which here stands for f x f, could be 
e 


confused with fe f.) 
In the alternative notation we would have 


b b 
[ arenyax or n| y? dx 


where y = f(x). 

There are two similar exercises to follow in which you are asked to find 
volumes of solids of revolution. If you are satisfied that you can cope with 
this type of exercise already, confirm this by tackling just one of the two 
and then moving on to section 9.3.2. 


Exercise I Exercise 1 
(3 minutes) 


Find the yolume ofa bucket of circular cross-section with the dimensions 
shown. | 
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Exercise 2 


Find the volume within the surface formed by rotating the graph of 
xe— x? (xe [1, 2]) 


about the x-axis. i} 


9.3.2 Averages 


In section 9.1 we looked at an example of an excavation problem, and 
determined the volume of earth to be removed. We did this by determining 
a cross-sectional area using given measurements. In actual practice, 
when viewing a site for building operations, a builder will often mentally 
“add a bit here” and “take a bit off there” and take an average depth to 
estimate the cross-sectional area. For example, faced with the site 
illustrated below a builder might assume that the area of the red region 
was the same as the area of the black region, and that the cross-sectional 
area required was the same as the area of the rectangle OABC. 


]rea 


Ground surface 


. 


He would probably call the height of this rectangle “the average height” 
of the site. If we suppose that the curve is the graph ofa function f between 
aand b relative to an appropriate axis, then the cross-sectional area is, 
of course, 7 


[is 
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Exercise 2 
(2 minutes) 


93.2 


Main Text 


(continued on page 45) 


Solution 9.3.1.1 


» 


The equation of the bounding curve is given by 
y-10  x-0 
15-10 25-0 


x 
y=5t 10 


[See Refresher Booklet 3 if you cannot follow this.) Therefore the volume 
expressed as 


rE) x 2 
a wo (a o} 
° 5 


1 (253-03) (25? — 0?) 
{is — 100(25 — of 


19 x 625 
EB) 3 


12 440 cm3 


= 12.44 litres (to 4 significant figures) a 


Solution 9.3.1.2 


2 
Volume = nf xtr— (x?)? 
4 


_ n(25 — 15) 


Solution 9.3.1.1 


Solution 9.3.1.2 


FM 9.3.2/9.3.3 


The base of our rectangle is then of length b ~ a, and if its height is H, 


(continued from page 43) 
we have 


ah 
average of f(x) over [ab] = H = rs | 7 if Definition 1 
= dd, * 
which we can abbreviate to: 
i b 
average = 
se= 7 i f 
Exercise | Exercise 1 
Find the average of f(x) over (0, 4], where f is the function: me) 
xt— x? (x € [0, 4]) 
| 
9.3.3 Velocity and Distance 9.3.3 


A hiker goes on a walk lasting four hours, excluding rests. From past Main Text 
experience he reckons that he walks at an average speed of 5 km/h. You aaa 
conclude that his walk was approximately 20km long. This was all 
straightforward since 


speed 


tres 


time 


speed x time = distance 


you assumed that the speed was constant, and so the calculation was 
simple. By plotting speed against time on a graph, we see that the area 
beneath the graph (the area of a rectangle in this case) represents the 
distance covered. 


Similarly, if the speed is not constant, the distance covered is also repre- 
sented by the area under an appropriate graph. We can see that this 
is so by dividing the time interval into n equal sub-intervals and then 
assuming that the speed is constant over each of these time-sub-intervals. 
The distance travelled during each sub-interval is then represented by 
the area of a rectangle, and the total distance covered is the sum of the 
areas of the n rectangles; this will be a close approximation to the area 
under the graph when n is large. 

Here we prefer to use “velocity” rather than “speed” because velocity 
means speed in a known direction. When we consider motion ina directed 
straight line, we take velocity to be positive if its direction is the same 
as that of the line, and negative if it has the opposite direction. 

The calculation of the distance travelled is not so easy as it was for the 
case of the hiker because the required area is given by a definite integral. 
We must be careful of the physical interpretation of “distance” in relation (continued on page 46) 
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Solution 9.3.2.1 


4 1/43 — 0 
Average i x x? = | = 54 | 
4-0], 


(continued from page 45) 


to the definite integral when the value of the velocity function becomes 
negative. In any particular context does it mean the total distance travelled 
or the distance that the object is from its starting point? This difference 
of meaning is illustrated in the following example. 


Exercise J 


A ball, thrown vertically upwards with an initial velocity of 20 m/s, has 
a velocity at time ¢ seconds given approximately by 


v(t) = (20 — 102) m/s 


Determine 


(i) the number of metres the ball is above the ground after 3 seconds; 
(ii) the number of metres that it has travelled in that time. 


(Exercise 9.2.2.2, parts (iv) and (v), page 33, will give you a hint how to 
approach this problem if you are in difficulty.) a 


9.3.4 Some Other Applications 


So far in section 9.3 we have looked at the applications of definite integrals 
to volumes of solids of revolution, averages, velocities and distances. 
Two engineering applications are very briefly outlined below. If you are 
not familiar with them, pass them by; they are merely introduced for the 
benefit of those who understand their background. 


Calculation of the total pressure on any surface immersed in water, 
such asa lock-gate or a dam, requires an integration, over the appropriate 
area, of the pressure at each point. The pressure in water increases linearly 
with depth. Ifa dam were triangular, this would lead to a definite integral 
of the form 


feos 


where z is the variable representing the depth below surface. 


Solution 9.3.2.1 


Exercise 1 
(3 minutes) 


9.3.4 


Discussion 


To determine the output of, say, a steam engine requires the evaluation 
of a definite integral of the type 


Puen 


where v is the volume, p the pressure in the cylinder and f is a known 
function which maps » to p. 


There are many more examples of the use of the definite integral in science 
and engineering, and in your further studies you will doubtless come 
across some of them. It is also interesting to note that a certain electrical 
process can be represented by a definite integral and that this fact can 
be used to advantage in the reverse sense : that is, we can design an electrical 
instrument which will evaluate definite integrals. For the moment we 
will use the word integrate to mean evaluate a definite integral. It will 
have a wider meaning in Unit 13. Integration IJ, but this one will suffice 
for the present. There is also a mechanical instrument, called a plani- 
meter, with which we trace out the boundary of the required area and 
which then records this area. If you think about it, a rain-gauge is an 
integrator. It is even an integrator if you don’t think about it*! It integrates 
the rate of rainfall over a specified period. And so also for that matter 
is a barometer —it adds up, or integrates, the weights of all the air 
molecules above it, the density (= mass/volume) at each level in the 
atmosphere being represented by quite a complicated function. 

Finally, look up the definition of “integral” in a dictionary. It is one of 
the words which is given a special mathematical meaning and yet retains 
much of its everyday sense. 


9.4 MORE REFINED APPROXIMATION METHODS 
FOR DEFINITE INTEGRALS 


9.4.0 Introduction 


In this unit we have defined the definite integral of a function as the limit 
of a sequence, and we have shown that for certain cases (simple poly- 
nomial functions) we can evaluate this limit. In Unit 13, Integration I, 
we shall obtain a general result which will allow us to extend the set of 
functions for which this evaluation is practicable. But even then there will 
still remain many functions for which that general process is impracticable, 
for example, 


1 
x-——- ————, 
} ° Jeri 
or for which it is unnecessarily complicated, for example, 
0.5 x? 
— 


z 
0 Lx 


When we meet integrals such as these in practical work, we usually need 
a numerical answer correct to some given accuracy, so we can often 
usefully return to our original approximation processes (or variations of 
these). Thus we obtain an estimate of the limit to a given accuracy, rather 
than an exact formula for the limit. . 


* We add this comment to our colloquialism because our next unit is on logic! 
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Introduction 


Solution 9.3.3.1 


Distances are represented by arcas on, the following diagram, as in the 
previous diagram for the simple example of the hiker. 


f velocity v 


-104 


We could in this case use the formula for the area of a triangle, but we 
will apply integration techniques to illustrate the general method. 

The graph crosses the t-axis at t = 2, indicating that after 2 seconds the 
velocity of the ball is zero. Subsequently the velocity is negative, i.e. the 
ball is returning to earth, so the ball will be at its highest point after 
2 seconds. To calculate the distance travelled up to this time we evaluate 


2 
i t-— (20 — 10r) = 20 
° 


So the ball is 20 m above its starting point after 2 seconds. 


1 


5m 


O@! 


20m| 


In the next second it travels a distance given by 


3 
i t-— (20 — 10t) = —5 
2 

the negative sign indicating, as expected, that the ball is returning towards 
the ground during this period. 


The answers are, therefore: 

(i) (20 — 5)m = 15m, 

(ii) (20 + 5)m = 25m. 

The answer to (i), 15 m, is the result of the evaluation of 


3 
f (¢-— (20 ~ 102) 
0 


that is, the definite integral represents the distance from the starting 
point. Notice the more general point here that, if we form the definite 
integral over a domain for which the image changes sign, we need not 
split the domain up into sub-domains in which the numerical value is 
wholly positive and wholly negative. We need only do this in cases where 
the physical requirements of the problem demand it, such as distance 
travelled in part (ii) of the example above, or as in part (v) of Exercise 
O22. | 
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Solution 9.3.3.1 


9.4.1 The Trapezoidal Rule 


We return now to a copy of the figure we used in the example at the 
beginning of section 9.1.2, when we were finding the area of part of a 
parabola. You will remember that we found a good estimate of the area 
by using 


H(A, + a,) 


where A, = the sum of the n areas of the larger rectangles, and a, = the 
sum of the n areas of the smaller rectangles where n = the number of 
intervals. 


treo 


What does this mean geometrically? Half the sum of one larger rectangle 
and one smaller rectangle (an example is outlined in the figure) equals 
the area of the trapezium ABCD (with upper boundary, CD, shown by 
the dashed line in the diagram). Thus the approximation we made earlier 
was equivalent to drawing the set of dashed lines (one for each interval) 
as the upper boundary to the area. That is, we approximate to the total 
area by the sum of the areas of the trapezia thus constructed. 


We now turn to the case where f(x) is non-negative in [a, b]. 


b 
Suppose we wish to find i J where the graph of f is given below: 
a 


f(x) 


»* 


We construct the set of dashed lines just as in the case of the parabola. 
Suppose the ordinates of points on the graph at 


a,a+h,a+2h,....a+nh=b 
are 


Yor ¥irYar-++>¥n Tespectively . 


where 
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9.4.1 


Main Text 


FM 9.4.1 


We have, from Exercise 9.1.0.1 on page 3: 
(area of first trapezium on left) = Ho t+ 1) xh 
(area of second trapezium from left) = Hy, + y2) x h 
and so on, until : 
(area of last trapezium) =Hyn-1 + Yn) XA 
Therefore, adding these equations together, we get: 


the total area of all the trapezia 
I 
= 50 + 21 + ya +o + Wana + Yad 
and therefore 
ee 
f SL = 5Vo + 291 + Wo +9 + na + Ya) 


where = means “‘is approximately equal to”’. 

This is known as the trapezoidal rule for evaluating definite integrals.* Definition 1 
Asindicated in Unit 2, Errors and Accuracy, we must now ask the question : ore 
“How approximate is ‘approximate’?” 

To answer this question we again consider the parabola. We use the 

fact that, for the rectangle approximation, provided the graph of the 

function slopes downwards over the whole interval or upwards over the 

whole interval, the difference, A, — a,, is the area of the largest rectangle. 

(See page 16.) In the case of the parabola: 


> 
| 
Q 


An — , = hi f(b) — f(a)| where h = 


The absolute error bound of the best estimate of the area, (A, + a,), is 
(A, — a,) which is 


4 f(0) — f(a| in the case of the parabola. 


The thing to notice about this formula is that, since f(b) and f(a) do 
not depend on the number of rectangles we choose to take in our approxi- 
mation, the absolute error bound is proportional to h, the width of each 
sub-interval. 

We return to the general case. The difference between this case and that 
of the parabola is that the graph of the function fno longer slopes down- 
wards over the whole interval or upwards over the whole interval. (See 
the second figure on page 49.) However, we have only to split [a,b] 
into smaller intervals in each of which the slope has the same sign 
throughout, and to apply the above argument to each of the smaller 
intervals. Since the “rectangle” method is equivalent to our trapezoidal 


* See note on page 54. 
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rule, this means that the best we can say at this stage about the absolute 
error bound for the trapezoidal rule is that it too is proportional to h. 


In fact, it is usually more accurate than the rectangle method. 


Example | 


How many intervals would you need to take, at most, to evaluate 


: 1 
x— 
4 x4 


to an accuracy of two decimal places using the trapezoidal rule? 


; ‘ 1 ‘ 
In practice, when evaluating Sar ve would round off the image 


+0 
values to some convenient number of decimal places, so that the error 
bound on each image will be ¢, say. To how many decimal places must 
the images be calculated to ensure that the integral is accurate to two 
decimal places, and how many intervals will then be needed to ensure 
that the overall accuracy is to two decimal places? Li 


Solution of Example 1 


An absolute error bound using the trapezoidal rule is 
h hy h 
sit) — fl =SB= 3 

where 


fix— 


b=landa=0 


x41 
(f(x) decreases as x increases in [0, 1].) 


For accuracy to 2 decimal places the absolute error bound must be less 
than or equal to 0.005 = 5 x 1073. 


Therefore 
h 
=<5x 10% 
ri x 
so 
h<2~x 107? 


and then, if the number of intervals is n, we have 


- 1 
na2= 425350 


Therefore a minimum of 50 intervals is required to guarantee the required 
accuracy (with our present knowledge about the accuracy of the 
trapezoidal rule). 


We have 
2 hy 

[ s=300 + 2yy +00 + Waa + Ya} 
On the right-hand side there are 2n values of f(x) inside the brackets. 
If each ordinate has an inherent error e, the total error on the right-hand 
side will be (see Unit 2, Errors and Accuracy) 

h 

= x 2ne = nhe = (b — ade 

2 
that is, a total error of ¢ in this particular case because (b — a) = | (note 


that this is independent of the number of intervals). 


SI 
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Example | 


To be really sure that our total error from the use of the trapezoidal 
rule and the inexact data does not exceed 0.005, we can use 100 intervals 
{the absolute error introduced by the trapezoidal rule from this is then 
< 0.0025) together with data accurate to 3 decimal places (error from 
this is then < 0.0005) to give a total possible error of 0.0030; but in fact 
50 intervals and data accurate to 3 decimal places would almost certainly 
suffice. ‘i | 


Exercise | 


Repeat the last example with the definite integral 
2 
i) x-— exp (— x?) 
0 


with a required accuracy of 3 decimal places. [e~* = 0.0183]. a 


9.4.2 Simpson’s Rule 


The trapezoidal rule gave an approximation to the area under a curve by 
using a set of straight line boundaries. One obvious way to improve 
the quality of the approximation is to take some account of the curvature 
of the boundary. In deriving the trapezoidal rule we took each pair of 
consecutive points (a + mh, y,,) and (a + (m+ 1)h, Yna1),m=0,...,2—1, 
found the straight line which passed through them (although we did not 
specify its equation since we already knew the area of a trapezium) and 
used this as the upper boundary of the area. We now introduce (no more 
than that) Simpson’s rule for the area under a curve. 


In Unit 4, Finite Differences, when we were studying polynomial inter- 
polation (section 4.3.2), we found that we could often interpolate (ie. 
find values between two tabular points) with increased accuracy if we 
increased the degree of the interpolating polynomial. We have an equi- 
valent situation here. The linear interpolation rule uses the same degree 
-of approximating polynomial, that is 1, as the trapezoidal rule (in which 
we replace the curved boundary by a set of line segments). The quadratic 
interpolation rule uses the same degree of approximating polynomial, 
that is 2, as Simpson’s rule (in which we replace the curved boundary by 
parabolic segments). The basic element of area is now the one shown in 
the diagram covering two intervals with a parabolic upper bounding 
surface. Thus we now need to split the total interval into “2-interval” 
or ‘‘3-point” subsets. 


What implication does this have on the number of sub-intervals we use? 
The answer to this is that the number of sub-intervals we use must now be 
even. 
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Exercise 1 
(5 minutes) 


9.4.2 


Main Text 


The bare outline of the basic steps in the argument are continued in the 
following text. You may like to check some of the steps and fill in the 
detail if you have the time; otherwise, read this section through quite 
quickly. 


Consider a “‘3-point” subset in which the 3 points to be fitted on the given 


curve have co-ordinates (~—h, yo), (0, y,), (h,y2). The approximating 
quadratic polynomial function has the form 


fix ax* +ayx+a  (xe[—A,h)) 
It can be shown: 
A 
(i) by evaluating | J, that the area beneath the graph of fis 
-h 


2a,h? 


+ 2agh 


(ii) by solving three simultaneous equations (or by using Lagrange’s 
interpolation formula: see Unit 4, Finite Differences) that 


a=) 
_ ¥2— Yo 
a Oh 
a, = 207 2H t+ V2 
. 2h? 


(iii) by substituting the appropriate results in (ii) into the result in (i) 
that the area beneath the graph of fis 


h 
zo + 4y1 + yo} 


This result gives the approximation to the area beneath the original 
given curve in any “3-point” interval in which the ordinates are yo, y; 
and y,. Thus in the next ‘‘3-point” interval with ordinates yz, y3, Ys 
(the interval in the figure shown dotted in black and red), the area approxi- 
mation is 


h 
ive + 4y3 + ya} 


froo 


ee 


Continuing in this way, we obtain the approximation (in the case of the 
function illustrated) to the total area as 


i 
give + 4ys + ya + Avy + 2y4 + dys + Yo) 
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(continued on page 54) 


Solution 9.4.1.1 Solution 9.4.1.1 
Your answer should contain the following: 
e~* decreases as x increases in [0, 2] 


2. 1-e%4 
n=— 


i 2 sx i077 = 2000 intervals 


If the error in each ordinate is ¢, then the error in the result is 2¢.- 


So we could use 2500 intervals (producing a possible error of 4 x 10~*) 
with data accurate to 4 decimal places (producing a possible error of 
Io 10"), a 


(continued from page 53) 


In fact, the general formula with n intervals (n even) is 


4 h 
[r= 500+ 4n + 2yp tore + 4Yn-1 + Yn) 
a 


This formula is known as Simpson’s rule. Definition 1 


It is generally true that you can obtain a better approximation (with the 
same number of intervals) using Simpson’s rule than with the trapezoidal 
rule, but again we need calculus to demonstrate this more explicitly. 
You will probably think that if we replace the boundary in each interval 
by fitting a cubic to four points we will get a better approximation than 
with the quadratic fitted to three points and using Simpson’s rule. This 
is not true, (although if we replace the boundary in each interval by 
fitting a quartic to five points, the accuracy does improve again). The 
trouble is that the complexity of the expression increases and one is in 
the usual ‘‘swings and roundabouts” situation. Only in cases when we 
get rather simpler formulation, such as when we use Gauss’s eleven 
point formula in which a tenth degree curve is used to approximate to 
the boundary in each interval, does the corresponding formula get a 
name, and even then it tends not to be used because the modern computer 
can perform the simpler repetitive numerical procedure of the trapezoidal 
and Simpson’s rule to the accuracy desired, by using many sub-intervals. 


Note 

We have obtained Simpson’s rule and the trapezoidal rule for {° f by con- 
sidering the area beneath the graph of f between a and b. For simplicity, 
we have considered the special case in which f(x) > 0 for all x in (a, 5). 
In fact, these rules apply when f(x) is not always positive in [a, b]. This 
can be seen by again considering areas (taking care of the signs) and 
slightly modifying our derivations. 
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9.5 CONCLUSION 


In this text we have used the historical link between area and integration 

in order to develop the definition of a definite integral as the limit of a 

particular sequence of sums. We have seen some of the wider applications 

of the definite integral (to volumes of revolution and averages, for example) 

and have found specific formulas for the definite integrals of some simple 

polynomial functions. If the function to be integrated is not a simple 

polynomial function, the definite integral (if it exists) can still be evaluated 

by 

either 

(i) using the more sophisticated specific formulas which will be derived 
in Unit 13, Integration H, if these are appropriate 

or 

(ii) finding an approximate value by the powerful tools introduced in 
the last section of this text. 


Acknowledgements 


Grateful acknowledgement is made to the Cement and Concrete Asso- 
ciation for permission to use the drawing of the Dollan Baths which 
appears on page 13, and also to the Mansell Collection for the portrait 
of Archimedes which appears on page 2. 


5S 


FM 9.5 


9.5 


Conclusion 
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